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Abstract 

We compute the production cross section of a top-antitop pair in association with a jet at hadron colliders at next-to- 
leading order accuracy matched with parton shower algorithms to make predictions at the hadron level. The parton 
shower allows for including the decay of the top quarks at the leading order accuracy. We use a framework based 
on three well established numerical codes, the POWHEG BOX, used for the calculation of the cross section, HELAC, 
which generates the matrix elements for the Born-level, real emission and the virtual part, and finally a parton shower 
program, such as PYTHIA or HERWIG, which generate the parton-shower and hadronization. 
PACS: 12.38.-t, 13.87.-a, 14.65.Ha 
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1. Introduction 

With the startup of the LHC, high energy particle 
physics entered a new era. At higher energies, measure- 
ments with higher precision become available, which 
poses new demands to the theoretical predictions: the 
corresponding cross sections are needed beyond lead- 
ing order (LO) accuracy even for large multiplicity final 
states. By now standard techniques exist mill for com- 
puting the next-to-leading order (NLO) corrections to 
many phenomenologically interesting processes involv- 
ing four, or more hard objects (heavy particle or hard 
jet) in the final state lO] |4] |5] 0! . Despite the improved 
accuracy obtained by computing the cross sections at 
NLO, there is still a large gap between fixed order the- 
oretical predictions and data collected by the detectors. 
At fixed order we calculate only hard parton-level pro- 
cesses, while in experiments we observe hadrons. The 
common practice to fill this gap is the use of parton 
shower programs [X ^| which also include hadroniza- 
tion models. The advantage of these programs is the 
generation of unweighted events, which can be utilized 
for performing the same analysis as on the collected 
data, allowing for a direct comparison of theory and 
experiment, or predicting the Standard Model back- 
ground. However, these programs catch only the impor- 
tant features of small angle radiation off" partons, and the 



distributions of observable quantities are not expected 
to give a good description in the regions dominated by 
large-angle hard emissions. 

Until recent years, these two main approaches were 
used separately for making predictions. Merging NLO 
computations with parton showers was pioneered by 
the MCatNLO project ||9). By now all interesting 
2 — > 2 processes are included in the MCONLD code fTOl. 
Another method for merging NLO computations with 
parton showers, which produces only positive weight 
events, was developed in refs. iTT T?]. The latter pro- 
cedure was later implemented in the POWHEG BOX [ISQ 
The POWHEG BOX can almost be considered a black box 
that requires matrix elements as input and produces un- 
weighted events in the form of Les-Houches accord files 
lfT4l as output. These events can be processed with the 
POWHEG BOX for generating the showered events for fur- 
ther analysis. 

In this letter we show the first application to a 2 — > 3 
process of the combination of the POWHEG BOX and the 
HELAC [2] frameworks for producing showered events 
of the 1 1 -H jet final state that can be used to make distri- 
butions with correct perturbative expansion up to NLO 
accuracy. Due to the large collision energy at the LHC, 
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Table I: Flavour stmctures of the Born processes, ^ = u, d, c, s, b. 
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Table 2: Flavour structures of the real-emission processes, q, q' = u, 
d, c. s, b. 

tt pairs with large transverse momentum will be copi- 
ously produced and the probability for the top quarks to 
radiate gluons will be sufficiently large to make the tt 
+ jet final state measurable with high statistics. There- 
fore, we make first predictions for such events at the 
TeVatron and the LHC. A more detailed analysis will 
be presented elsewhere. 

2. Method 

The cornerstone of our program is the POWHEG BOX 
ifTSll framework, that uses the FKS subtraction scheme 
ifTSll for the NLO calculation. The POWHEG BOX re- 
quires the following input: 

• We use flavour structures given in Tables [T] and |2] 

• We generate a Bom phase space of a massless and 
two massive momenta using two two-particle in- 
variants and three angles. 

• We use HELAC-Dipoles |]T6l to calculate all the 
tree-level helicity amplitudes for the Born subpro- 
cesses tlggg — > and tlqqg — > and the real 
emission subprocesses tlgggg — > 0, tlqqgg 

and ttqqq' q' — > 0. (We define the correspond- 
ing crossing symmetric amplitudes for all incom- 
ing momenta and cross into the relevant physical 
channels.) 

• For the colour-correlated squared matrix elements 
of the Born subprocesses we use HELAC-Dipoles. 

• We use the polarization vectors to project the he- 
licity amplitudes to Lorentz basis for writing the 
spin-correlated squared matrix elements. 

• Finally, we obtain the one-loop corrections to the 
Born subprocesses utilizing the HELAC-Oneloop 



implementation |l2] (TT] [HI of unitary-based nu- 
merical evaluation of one-loop amplitudes lfT9ll20l 

ED EH Ha El |25] US. 

Wifli this input POWHEG BOX can be used to generate 
hadronic events. One may choose any parton shower 
(PS) Monte Carlo program for generating parton show- 
ers, decays of heavy quarks and hadronization. There 
is one important point in choosing the PS. We gener- 
ate events with hardest emission measured by the trans- 
verse momentum of the emission. If the ordering vari- 
able in the shower is different from the transverse mo- 
mentum of the parton splitting (for instance, the angular 
ordered showers in HERWIG), then the hardest emission 
is not necessarily the first one. In such cases the HERWIG 
discards shower evolutions (vetoed shower) with larger 
transverse momentum in a subsequent splitting than that 
in the real emission correction. In addition, a truncated 
shower simulating wide-angle soft emission before the 
first emission is also needed in principle, but its effect 
was found small ll27l . As there is no implementation of 
truncated shower in HERWIG using external LHE event 
files, the effect of the truncated showers is absent from 
our predictions. 

3. Checks 

In order to ensure the correctness of the calculations 
we performed the following checks relevant to any fixed 
order calculation at the NLO accuracy: (i) Compared 
the cross section at LO to the prediction of the public 
code MADGRAPH ||29l and found complete agreement, 
(ii) Checked the virtual correction obtained from the 
HELAC-Oneloop program in several randomly chosen 
phase space points to that obtained from the implemen- 
tation in the PowHel (= PDWHEG+HELAC) program, (iii) 
Checked in several randomly chosen phase-space re- 
gions that the ratio of the soft- and collinear limits of 
the real-emission matrix elements and subtractions tend 
to one in all possible unresolved limits. 

There is an important technical issue related to the 
way of calculation organized in the POWHEG BOX. The 
selection cuts are applied on the events obtained after 
hadronization. However, when computing the 1 1 + jet 
production cross section at fixed order, the cuts are ap- 
plied at the parton level. At LO this means a cut on 
the transverse momentum of the only massless parton in 
the final state. At NLO the virtual contribution has the 
same event configuration as the Born one, but the real 
emission contribution has two massless partons in the 
final state, that have to be combined into a jet before the 
physical cut can be applied. In the POWHEG BOX such 
a separation of the real and virtual contributions is not 
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Table 3: Dependence of the NLO cross section on the technical cut 



possible because the event-generation starts with an un- 
derlying Born configuration from which further parton 
emissions are generated. In order to make the parton- 
level calculation finite, we can apply a technical cut on 
the transverse momentum of the single massless par- 
ton in the Born confiuration. With a given set of se- 
lection cuts, one has to check that the chosen technical 
cut is sufficiently loose such that it does not influence 
the physical cross section. Typically we find that for 
jet transverse momentum cuts of several tens of GeV, 
a several GeV technical cut on the transverse momen- 
tum of the massless parton at Bom level is sufficiently 
loose. Another way of treating the same problem, also 
implemented in the POWHEG BOX, is to use a suppres- 
sion factor on the underlying Born configuration ll28l . 

4. Comparison to predictions at NLO 

The first calculation of the tt + jet production 
cross section was computed by Dittmaier, Uwer and 
Weinzierl [30. .31.1 . In order to further check our pro- 
gram, we computed the production cross section at 
NLO accuracy using the same physical parameters as in 
Ref. [3 LI- Due to the technical cut mentioned in the pre- 
vious section, the PowHel framework is not optimal for 
a fixed-order computation, nevertheless our prediction, 
^NLO - (1 78 + 0.01) pb is in agreement with the cross 
section quoted in Ref. ED, o-Nlo = (i 791 + 0.001)pb, 
within the uncertainty of our integration. Our prediction 
is independent of the technical cut below p^f- < 5 GeV 
as shown in Table |3] 

In order to check the predictions obtained with Born- 
suppression, we computed the distributions published 
in Ref. OTll at NLO accuracy and we found agreement. 
Examples are shown in Fig. [T] for the case of the trans- 
verse momentum and rapidity distributions of the jet. 
The lower panels show the ratio of the PowHel-NLO 
predictions to the predictions of Ref. 1 3 1 1 . The error 
bars in the lower panel represent the combined statis- 
tical uncertainty of the two computations. 

We also compared distributions obtained from LHE 
events, including the first radiation only, to predictions 
at NLO. For the distributions of the transverse momenta 
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Figure 1: Transverse momentum and rapidity distributions of the jet. 

of the jet (see also in Fig.[T]) and the top as well as for 
the rapidity distribution of the top we found agreement. 
The rapidity distribution of the jet is more central from 
the LHE events than from NLO. 

5. Effects of decays and shower 

The production of 1 1 + jet final state at the NLO ac- 
curacy together with decay of the heavy quarks in the 
narrow-width approximation (at LO accuracy) has been 
published recently by Melnikov and Schulze in |32l . In 
our NLO+PS computation decays of heavy quarks are 
implemented in the PS, therefore, spin correlations are 
not included. In contrast, the narrow-width approxima- 
tion allows for taking into account the spin correlations. 
Thus, in order to see the effect of the parton shower, we 
first generated distributions without the shower, but with 
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decays (we just included on-shell decays of t-quarks, 
and further decays of their decay products, if unstable, 
turning off any shower and hadronization eifect, marked 
as 'Decay'), then with the full shower Monte Carlo 
(marked with the name of the SMC). We compared 
the total cross section as well as several distributions 
to those predictions made for collisions at the Teva- 
tron, = 1 .96 TeV, valid at the NLO accuracy. We 
generated two million events with PowHel, which were 
showered with PYTHIA-6.4.25 ||33] and HERWIG-6.5.20 
0411 subsequently. For the comparison, we used the 
semileptonic decay channel and the following parame- 
ters and selection cuts from Ref. 1321 : (i) mass of the top 
quark nii - 172 GeV; all other Standard Model param- 
eters as implemented in the PS programs, (ii) CTEQ6M 
parton distribution functions, (iii) A;^ -clustering algo- 
rithm with R - 0.5 and four-momentum recombination 
scheme |36|, (iv) /iR = /ip = nit, (v) > 20 GeV, (vi) 
> 20 GeV, (vii) p{ > 20 GeV, (viii) \yj\ < 2, (ix) 
minimum five jets, and (x) Hj^ > 220 GeV, where //^ 
is the scalar sum of transverse momenta in the event, 
= + Ef'' + Zjp{ ■ In addition, if the final 
state after these selection cuts contained one or more 
charged leptons, we rejected the event if the transverse 
momentum of this lepton was above 20 GeV. This latter 
requirement is not needed in a fixed order calculation, 
but necessary in ours to select the semileptonic channel. 
The technical cut was chosen to p^f = 5 GeV. 

The predicted SMC cross sections are very sensitive 
to the details of the analysis. We kept the leptons and 
neutral pions stable, while all other particles were al- 
lowed to be stable or to decay according to the default 
implementation in each SMC. Quark masses, as well 
as W, Z masses and total decay widths, were tuned to 
the same values in PYTHIA and HERWIG. On the other 
hand, each of the two codes was allowed to compute 
autonomously partial branching fractions in diff'erent 
decay channels for all unstable particles and hadrons. 
Multiparticle interaction effects were neglected (default 
in HERWIG). Additionally, the intrinsic pr spreading of 
valence partons in incoming hadrons in HERWIG was as- 
sumed to be 2.5 GeV. 

Considering this setup, we always found agreement 
between PYTHIA and HERWIG predictions within 3 %, 
which is also the effect of including versus neglect- 
ing negative weight events in the analysis. For in- 
stance, using our selection cuts and taking into account 
the negative weight events, we obtained the cross sec- 
tions o-P°"Hei+HERwiG ^ 146.9 ft and o-p°""^i+p™" = 
143.2 fb, while without the negative weight events, we 
obtain crP°"«^i+PYTHiA ^ ^^-j corresponding 



value for the PowHel-i-decay case is crP™H'='+'i'='^^y - 
144.2 fb (with negative weight events included). These 
numbers cannot be compared directly to the fixed-order 
prediction o-Nlo ^ 33.6 fb quoted in Ref. (321 for two 
reasons. On the one hand in Ref. ['321 only one lep- 
ton family was considered in the decay of the t-quarks, 
while our prediction contains all three families. We 
checked that taking into account only one lepton family 
in the decay we obtain a factor of three reduction of the 
cross section as expected. On the other hand the authors 
of Ref. 1 32 1 also observed that there is a large contri- 
bution to the cross section from the emission of a hard 
jet from the top decay products (estimated an additional 
60 % at LO |331), which is included in our calculation, 
but not in their value. As this effect is not known at the 
NLO accuracy, in order to compare only the shapes of 
distributions with only decays included, we multiply the 

NLO predictions with r = crPowHel+decay/^NLO ^ 4 29 

(shown as 'NLO-(-decay' in Figs. |2] and |3]l. The lower 
panels show the ratio of the various predictions to the 
PowHelH-PYTHIA one. In order to exhibit the size of 
the statistical uncertainty (corresponding to two mil- 
lion LHE events), avoiding at the same time a very 
confusing plot, we show the uncertainty of only the 
PowHelH-decay prediction with errorbars. 

In Fig.|2]we compare the transverse momentum and 
rapidity distributions of the antilepton at several dif- 
ferent levels. We observe on these plots some gen- 
eral features: (i) the two PowHel-t-SMC predictions are 
very close except in bins with low statistics; (ii) the 
PowHel+decay predictions are very close to the NLO 
ones in the central rapidity region and for the whole 
p_i_ range. Looking more closely, we find that the spin 
correlations make the NLO rapidity distribution slightly 
wider. The addition of the parton shower makes the ra- 
pidity distribution a little even more central due to soft 
leptons emitted by the shower in central regions. (For jet 
rapidities, not shown here, the NLO and PowHelH-decay 
predictions coincide, but the shower effect is much more 
pronounced.) The p^-distributions of the leptons be- 
comes much softer for the same reason. The same ap- 
plies to the spectra of the jets. 

We find even larger shower effects in the comparison 
of the i/j^ -distributions in Fig.|3]at the decay and SMC 
levels. The shower makes the distribution softer, readily 
understood as the effect of unclustered soft hadrons in 
the event, that appear only in the shower. 

6. Predictions for the LHC 

We now turn our attention to the LHC and make some 
predictions for the inclusive tt -1- jet production at the 
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Figure 2: Transverse momentum and rapidity distributions of the an- 
tilepton. 



low-energy run, yfs -7 TeV in the dileptonic final state 
channel. We apply the following selection criteria: (i) 
at least three jets are reconstructed with the anti-A;j_- 
clustering algorithm with R - 0.5 and four-momentum 
recombination scheme |37|, (ii) > 30 GeV, (iii) 
ly^l < 2.5, (iv) Ef'' > 30 GeV for e+e" and fi+fi' 
pairs, while /J™'**^ > 20 GeV for e'^fi^ pairs, (v) , 
> 20 GeV for exactly one and one 
For default scales we used two different choices: (i) 
the mass of the t-quark, nit, and (ii) the transverse 
mass of the harder top, yUR = /ip = m^, where mj_ - 

^ \ We expect the latter scale better 
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interpolates between near-threshold and hard events. 

In Fig.|4]we plot the transverse momentum distribu- 
tions of the hardest, second hardest and third hardest 
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Figure 3: Distribution of the scalar' sum of transverse momenta. 

jet. These /7x spectra are insensitive to the version of 
the parton shower within the statistical uncertainty of 
the computations, which shows that the effect of the 
missing truncated shower must be small. Also they are 
rather robust against the choice of the default scale (2- 
6 % variation, not shown here), suggesting small scale 
dependence in general, but we shall study that in a sepa- 
rate publication. The same features are also true for the 
rapidity distribution of the antilepton as seen in Fig.|5] 
where we also exhibited the prediction at the decay 
level. The lower panel shows the ratios of the predic- 
tions to the PowHel+HERWIG case. The error bars rep- 
resent the statistical uncertainty of the latter only. We 
find large (almost 20 %) and almost uniform effect of 
the shower and hadronization. In the case of the trans- 
verse momentum distribution of the antilepton the var- 
ious predictions agree over the whole spectrum except 
that we see a large increase from the decay level to the 
full SMC at small p^, see Fig.|6] We attribute this in- 
crease to the numerous secondary antileptons generated 
in the hadronization phases. 

Finally, we plot the invariant mass distribution of the 
i^l'^ pairs in Fig.|7] Here again the full SMC predic- 
tions are all the same. During hadronization additional 
(anti)leptons with > 20 GeV may appear and such 
events are dropped due to our selection cut (v), result- 
ing in a softer spectrum. 

7. Conclusions 

In this paper we interfaced the POWHEG BOX with 
the HELAC framework to perform NLO calculations 
matched with parton showers and hadronization in a 
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Figure 4: Transverse momentum distributions of the first, second and 
third hardest jet. 

quite general and semi-automatic way. The latter means 
that the necessary ingredients for the POWHEG BOX can 
be taken from the HELAC framework without any further 
computations. We presented the feasibility with a non- 
trivial process, namely inclusive 1 1 + jet production and 
decay, and we found reliable results. We employed de- 
cays as implemented in standard PS Monte Carlo pro- 
grams. We leave the extension to decays included in the 
hard matrix elements for a future study. 

We emphasize that the necessary virtual emission was 
calculated by a general numerical method which can be 
used for further processes. Due to the general nature of 
our framework including further processes is feasible. 

Using the PowHel framework we produce several 
million unweighted events at the hadron level readily 
available for analysis. These events can be used to 
produce distributions that are correct at NLO accuracy 
when expanded in the strong coupling. Our analyses 
clearly show the importance of the full SMC. There are 
certain regions in the phase space, where even a NLO 
accuracy is insufficient. A singular example is the i/^ 
distribution which shows significant softening over the 
whole kinematic range. 

In preparing this letter we learnt about a similar work 
in progress by Alioli, Moch and Uwer, presented at the 
Heavy particles at the LHC workshop, Zurich, 20 11. 
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Figure 5: Rapdity distribution of the antilepton. The lower plot also 
includes the ratio of the cross section obtained with fi = to that 
obtained with /i = (PowHel+PYTHIA). 
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with /J = m, to that obtained with /i = m± (PowHel+PYTHIA). 
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